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Abstract 

I present a unified approach to calculating the next-to-next-to-leading order (NNLO) 
soft and virtual QCD corrections to cross sections for electroweak, Higgs, QCD, and 
SUSY processes. I derive master formulas that can be used for any of these processes 
in hadron-hadron and lepton-hadron collisions. The formulas are based on a unified 
threshold resummation formalism and can be applied to both total and differential cross 
sections for processes with either simple or complex color flows and for various factor- 
ization schemes and kinematics. As a test of the formalism, I rederive known NNLO re- 
sults for Drell-Yan and Higgs production, deep inelastic scattering, and W~^^ production, 
and I obtain expressions for several two-loop anomalous dimensions and other quantities 
needed in next-to-next-to-leading-logarithm (NNLL) resummations. I also present new 
results for the production of supersymmetric charged Higgs bosons; massive electroweak 
vector bosons; photons; heavy quarks in lepton-hadron and hadron-hadron collisions and 
in flavor-changing neutral current processes; jets; and squarks and gluinos. The NNLO 
soft and virtual corrections are often dominant, especially near threshold. Thus, a unified 
approach to these corrections is important in the search for new physics at present and 
future colliders. 



1 Introduction 



Progress in theoretical particle physics, from electroweak theory and Quantum Chromodynam- 
ics in the Standard Model to supersymmetry and beyond, often involves the comparison of 
predictions of theory with data from high-energy hadron-hadron and lepton-hadron collisions. 
Recent theoretical advances include an array of resummations and a few next-to-next-to- leading 
order (NNLO) calculations [HE]- These advances are necessary in many cases where next-to- 
leading order (NLO) calculations are not accurate enough, since higher-order corrections reduce 
the scale dependence and increase theoretical accuracy. 

Next-to-next-to-leading order calculations are technically very challenging and have been 
completed only for a few processes, including Drell-Yan [3 H] and Higgs production El EI 
and deep inelastic scattering [HI El- The corrections are usually split into hard, soft, and virtual 
parts, corresponding to contributions from energetic, soft, and virtual gluons, respectively. The 
soft and virtual corrections are an important component of the total result both theoretically 
and numerically. In fact, in some schemes and kinematical regions, e.g. threshold, they are 
the dominant part. Highlighting their importance is the fact that for Drell-Yan JU] and Higgs 
|llllT^ production the NNLO soft and virtual corrections were presented before the full NNLO 
result was calculated. 

We will see that there is a universality in the form of these corrections, which becomes 
more evident from the techniques of threshold resummations, which arise from factorization 
properties of the cross sections. Threshold corrections can be resummed to all orders, and 
finite-order expansions of resummed cross sections have provided us with many cross sections 
at NNLO and next-to-next-to- leading logarithm (NNLL) accuracy [13j. Thus it is a worthwhile 
aim to see if a unified approach can be given for the calculation of NNLO total and differential 
cross sections for any process in the Standard Model and beyond in hadron-hadron and hadron- 
lepton colliders and in various factorization schemes and kinematics. This will not only increase 
and deepen theoretical understanding but will also help avoid effort duplication in calculations 
of NNLO corrections for new processes. It is important to note that new particles, such as in 
supersymmetry, will likely be discovered near threshold where the soft and virtual corrections 
are important. The new unified approach to NNLO soft and virtual corrections is the topic of 
this paper. 

The calculation of cross sections in hadron-hadron or lepton-hadron collisions can be written 
schematically as 



where a is the physical cross section, (pj/hi is the distribution function for parton / carrying 
momentum fraction Xi of hadron /ij, at a factorization scale while is the renormalization 
scale. The parton-level hard scattering cross section is denoted by a, and s and are standard 
kinematical invariants. In a lepton-hadron collision we obviously have one parton distribution 
{i = 1) while in a hadron-hadron collision i = 1,2. We note here that a and a are not restricted 
to be total cross sections; they can represent any relevant differential cross section of interest. 

In general, a includes plus distributions Vi{xth) and delta functions S{xth) with respect to 
a kinematical variable Xth that measures distance from threshold, with / < 2n — 1 at nth order 
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in as beyond the leading order. These are the soft and virtual corrections. In single-particle 
inclusive (IPI) kinematics, Xth is usually called S4 (or S2), = s + t + u — "^^m?, and vanishes 
at threshold. Then 



S4 



S4 ' ' l+l VM2 



where A is a small parameter introduced in order to separate the hard Xth > A and soft Xth < A 
gluon regions, and is a hard scale relevant to the process at hand, for example the mass 
m of a heavy quark, the transverse momentum pt of a jet, etc. In pair-invariant-mass (PIM) 
kinematics, with the invariant mass squared of the produced pair, xth is usually called 1 — x 
or 1 — z, with z = Q'^/s — > 1 at threshold. Then 



ln'(l-^)l ln'(l-zl 



9{l-z-A) + - — -W+\A)6{1- z) . (1.3) 



The highest powers of these distributions at each order in a.^ are the leading logarithms (LL), 
the second highest are the next-to- leading logarithms (NLL), the third highest are the next- 
to-next-to- leading logarithms (NNLL), etc. These logarithms can be resummed to all orders 
in perturbation theory. By now there are several processes for which NLL resummations and 
NNLO-NNLL results (i.e. the NNLL terms at NNLO) have been presented fH^ . 

In this paper I will present master formulas for the NLO and NNLO soft and virtual 
corrections for any process in hadron-hadron and lepton-hadron collisions. In the next sec- 
tion, I present a threshold resummation formula, that builds on and unifies previous work 

I14[ [W\ IT7| I18j . I then present master formulas for the NLO and NNLO soft and virtual 
corrections in processes with simple color flow, that arise from the expansion of the resumma- 
tion formula and matching to NLO. The formulas cover both the MS and DIS schemes, and 
IPI and PIM kinematics. In Section 3, I present results for electroweak/Higgs processes as 
well as QCD and SUSY processes with simple color flows. I rederive known NNLO results for 
the Drell-Yan process and Higgs production, as well as for deep inelastic scattering and W'^'j 
production, thus obtaing expressions for two-loop anomalous dimensions and other quantities 
that are universal in quark-antiquark and gluon-gluon scattering and are needed for NNLL re- 
summations. I also present new results for the NNLO corrections for supersymmetric charged 
Higgs production, W, Z plus jet production, direct photon production, DIS heavy quark pro- 
duction, as well as single-top quark production mediated by flavor-changing neutral currents. 
In Section 4, I extend the master formulas to processes with complex color flows and I present 
results for heavy quark hadropro duct ion, jet production, and squark and gluino production. I 
close with a discussion of extensions to higher orders and conclusions in Section 5. 



2 NNLO master formula for soft and virtual corrections 
2.1 Soft corrections from threshold resummation 

In hadron-hadron collisions we study processes where partons fi collide and produce a speciflc 
flnal state. The partonic processes are of the form 

+ f2{p2)^F + X, (2.1) 
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where F represents a system in the final state, and X any additional allowed final-state particles. 
So F can represent a pair of heavy quarks, or a single heavy quark, or jet, or photon, a Higgs 
boson, squarks, etc. 

In lepton-hadron collisions the processes are of the form 



F + X . 



[2.2) 



In either case s = {pi + P2Y ^^'^ ^« usual Mandelstam invariants formed by the four- 

momenta of the particles in the scattering. 

Resummed cross sections have by now been studied for a variety of processes fTB]. The 
resummation of threshold logarithms is carried out in moment space. We define moments of 
the partonic cross section by (t{N) = J dz z^~^a{z) or by a{N) = j{dsi/s) e~^*'*/**o"(s4), with 
the moment variable. The resummed partonic cross section in moment space is then given 

by 



exp 



X exp 



exp 

,/2 



Tr [h {a.{^,l)) Pexp [ ^ T% {a,{^^'-) 



exp 



2d 



~S {as{s/Nf) 



X Pexp 



:2.3) 



The sums over i run over incoming partons: in hadron-hadron colisions we have two partons 
in the initial state, so i = 1, 2; in lepton-hadron collisions we only have one parton (and one 
lepton). The sum over j is relevant only if we have massless partons in the final state at lowest 
order. Clearly the second exponent is then absent for proceses such as Drell-Yan, Higgs, and 
top quark pair production. 

Equation ()2.3|) is actually valid for both IPI and PIM kinematics with appropriate definitions 
for Ni and Nj. In IPI kinematics Ni = N{—ti/M'^) for incoming partons i, and Nj = N{s/M'^) 
for outgoing massless partons j; here is any chosen hard scale relevant to the process at 
hand. The kinematical invariants ti are assigned through S^/ S = s^/ s — ~ Xi)ti/ s [T71 [12] 
with 5*4, 5*, the hadronic analogs of S4, s; note that the ti may include a —m? in case of massive 
particles. In PIM kinematics = Nj = N. Often the resummed cross sections are given with 
the choice = s in which case Nj = N even in IPI kinematics; here we keep our expressions 
more general. Also note that = Ne'^'^, with 7^ the Euler constant. 

The first exponent in Eq. ()2.3p resums the A^-dependence of incoming partons fl^ 1^ and 
is given in the MS scheme by [IHj 



dz- 




dfi 



{l-z)^s fl' 



[2A) 
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with = C/Ja,/7r + {as/'KfK/2] + ■ ■ -. Here C/^ = Cp = {N^ - l)/(2iV,) for an 

incoming quark or antiquark, and Cf- = Ca = Nc for an incoming gluon, with Nc the number 
of colors, while K = Ca (67/18 — 7r^/6) —5nf/9, where nf is the number of quark flavors. Also 
uj^ = {asMCf^ + {as/r^fvf + ■ ■ -. 

The second exponent resums the iV-dependence of any outgoing massless partons and is 
given by jTTl HH] 

E'^f^\N,) = / dz— / [as [Xs)) - B'^ - z)s)] - v, - zfs) 

JO i — Z [J(l—z)'' A '- 

(2.5) 

where B'j = {as/'n)B'f + [asl-nfB'f + ■■■ with B'f = 'iCpj^ and B'f = po/A. We will 
determine B'^^^ in Section 3. 

The 7^^"* in the third exponent are one-loop parton anomalous dimensions: 7^^-* = 3Cf/4: 
and 7^^-* = Pq/A for quarks and gluons, respectively; it is important to note that in this specific 
form that the resummed cross section has been written, with as the upper limit of the 
integral over dn''^ in E^^^\Ni), only the one-loop 7j appears in the third exponent. We also 
have defined 7'j/j as the moment-space anomalous dimension of the MS density (pi/i, minus it's 
one-loop and its A^-dependent two-loop components. This is again due to the specific form that 
we use for the resummed cross section. Thus, j -^^ is the A^-independent part of the two-loop 
anomalous dimension 7j/j [221 and is given for quarks and gluons by 

VSl ^ CI (1 - Ic. . fc) + C.C., (-?C3 . , 1!) ^ „,c. (-1 - , (2.6, 

and 

respectively, with ^2 = 7r^/6 and (3 = 1.2020569 ■ ■ ■. 

The (3 function in the fourth exponent is given by (3{as) = iidlng/d^i = — /5oas/(47r) — 
pia^JiAnf + ■ ■ -, with /?o = (HCa - 2nj)/3 and /3i = 34C1/3 - 2nf{CF + 5Ca/3). The 
constant da^ = 0, 1, 2 if the Born cross section is of order a^, a], a^, respectively. 

The trace appearing in the resummed expression is taken in the space of color exchanges. 
The symbols P and P denote path ordering in the same sense as the variable n' and against 
it, respectively. The evolution of the soft function from scale \fs/Nj to ^/s follows from its 
renormalization group properties and is given in terms of the soft anomalous dimension matrix 

[ini HEl Unj- In Eq- ()2.3j) we actually use r'5, which is given by Vs after dropping all 
gauge- dependent terms. We can do that because the gauge dependence has been shown to 
cancel out. At one loop, the gauge terms are of the form Cf^ ln(2i/"), where z/" = (fj ■ n)^/|np 
with Vi a velocity vector and n the axial gauge vector ^3 El . In processes with simple 
color flow, r5 is a trivial 1x1 matrix. For the determination of Ts in processes with complex 
color flow, an appropriate choice of color basis has to be made. For gluon-gluon scattering, the 
most complex color flow encountered, F^ is an 8 x 8 matrix 16 . The process-dependent soft 
anomalous dimension matrices, evaluated through the calculation of eikonal vertex corrections 
|15j . have by now been presented at one loop for practically all partonic processes ^H]- They 
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can be explicitly calculated for any process using the techniques and results in Refs. [131 IT^ . 
Work is currently being done on two-loop calculations of these anomalous dimensions [21] , but 
we note that we can extract the universal component of these anomalous dimensions for quark- 
antiquark and gluon-gluon initiated processes from the Drell-Yan and Higgs NNLO results as 
detailed in Section 3. In the color bases that we normally use, the soft matrices, S, are diagonal. 
At lowest order, the trace of the product of the hard and soft matrices reproduces the Born cross 
section for each partonic process. We also note that the matrices are in general not diagonal 
in the color bases that we use for complex color flows. If we perform a diagonalization so that 
the Ts matrices do become diagonal, then the path-ordered exponentials of matrices in the 
resummed expression reduce to simple exponentials; however, this diagonalization procedure is 
complicated in practice ^H] • A finite-order expansion bypasses the need for this diagonalization 
procedure. For the discussion below we expand T's as F'g = {as/7T)T'g''^ + (as/7r)^F'^'* + ■ ■ •. 
We will determine the universal component of ReF'^^'' — where "Re" stands for "real part 
of," for quark-antiquark and gluon-gluon processes in Section 3. 

Finally, we note that there are different ways of writing threshold resummation formulas 
that have been presented in the past for various processes, all consistent or equivalent at NLL 
accuracy; the expression presented here unifies those expressions for arbitrary processes and is 
superior in its simplicity and generality. 

In the DIS scheme, the resummed cross section may be written as 



where is the MS cross section in Eq. (jTHI), and % = (a,/7r)5j^^ + {as/T^fsf + ■ ■ -. 
For quarks, to which the scheme is usually applied, B^"* = ?)Cp/A. We will determine B^^ in 
Section 3. 

We first expand the resummed formulas in Eqs. (j2.3p and (j2.8j) for processes with simple 
color flow, i.e. when H, S and F5 are 1x1 matrices, to next-to-leading order and present NLO 
master formulas in the MS and DIS schemes, which reproduce the NLO results for a variety of 
processes. We then present NNLO master formulas for soft and virtual corrections in processes 
with simple color flow in both the MS and DIS schemes. 

2.2 NLO master formula for soft and virtual corrections - simple 
color flow 

2.2.1 MS scheme 

At next-to-leading order, the expansion of Eq. ()2.3p gives the NLO soft and virtual corrections 
in the MS scheme: 

^(1) ^ {c,V,{xth) + + ci 5{xth)} , (2.9) 

vr 

where is the Born term, 

C3 = E2C/,-E^/.' (2-10) 

i j 



(2 



6 



and 



C2 



Cf, + 2Cf, 6k In 



(2.11) 

where Rer'^^"* is the real part of the one-loop Ts after dropping all gauge-dependent terms, and 
5k is for PIM kinematics and 1 otherwise. We remind the reader that the sums over i run 
over incoming partons and the sums over j run over any massless partons in the final state 
at lowest order. For future use we will write C2 = C2 + T2, where C2 represents the scale term 
— J2i ln{ijLp/s) and T2 is the remainder. Also, ci — + 71, with 



E 



f. Ok In ( — ^ 



Cf.S 



-r. 



(1) 




(2.12) 



Ti denotes the terms in ci that do not involve the factorization and renormalization scales, 
and it can be read off from a full calculation of the NLO virtual corrections for any specified 
process; a formal expression for these terms is given in Section 4.1 in the more general case of 
complex color flow. 

We note that our formula passes a number of tests. As we will see, its predictions agree with 
exact NLO soft plus virtual results for all processes where those results are already available. 
Also the renormalization and factorization scale dependence in the physical cross section (after 
convoluting the partonic cross section with the parton distributions) cancels explicitly, i.e. 
da/d/ip = and da/d^R = at NLO. 



2.2.2 DIS scheme 

In the DIS scheme, the NLO soft and virtual corrections are 



'DIS 

with c'3 = C3 - EjC'/,, 



CniG = cr 



TT 



{cg Vi{xth) + C2 Vo{xth) + c[ S{xth)} , 



C,Sk^{^^B^ 



and 



Vm2 



(2.13) 

(2.14) 
(2.15) 



with C3, C2, Ci the MS results given in the previous subsection. For future use we will write 
C2 = -|- T2, and c'^ = c^'^ -|- T{ as we did for the MS corrections. Note that the changes 
in going from the MS scheme to the DIS scheme are all in the scale-independent parts of the 
Cj's, i.e. C2 and Cj* remain unchanged while T2 and Ti (and C3) are modified when changing 
schemes. The DIS scheme is normally applied to quarks. For quarks the term Df- in c[ is 
Dq — CfC,2 + (9/4)Cir. We note that our formula passes the same tests as we outlined in the 
previous subsection for the MS scheme. 
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2.3 NNLO master formula for soft and virtual corrections - simple 
color flow 

2.3.1 MS scheme 

At next-to- next-to-leading order, the expansion of Eq. ()2.3jl . with matching to the fuU NLO 
soft-plus-virtual result, gives the NNLO soft and virtual corrections in the MS scheme 



(2) 

' = a — a 



with 

,(2) 



a 



-Cg V^{Xth) + 



-C3 C2 - ^C3 + }^ Cf^ 



/3o 



(2.16) 



T^2{Xth) 



+ i C2 Ci - C2 C2 C3 + Cs C3 



|T, + §c.ln(^f^)+2Rerf-i:-f 



In 

2 V s 



'(2) , ,,(2) 



+ <! - I + 1(2^ + Cs C3 C2 - ^C4 + In ( ^ ) + 2 Rer'^^^ 5^ In ( — 



° In^ 



l--ln 



1 

2 

/3o 



1 

4^ 



3 
-( 

4 



2 

+ E 



Ti In 

(1) 



+ —6kT2 In 



4 



s 



+ 



16 



-^0 ( f^R 



+ /3i In 



/3o 



- 6k In 



E7l/1ln^)+EC^.^- 



6 VA/V 



+ 'T+ 2 



In^ 



(ill) - (ill) 



In 



In 

- In I — 

2 V s 



Vilf2 



+ In 



(ilfO 



+ E '^z?^ '^i^ In 

i 

+ Ei^Q.ln 



(ilfO 



(2.17) 
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where the last term R, for which a formal expression is given in Section 4.1 in the more general 
case of complex color flow, can only be known from a full two-loop calculation. We will derive 
in Section 3 the universal components of R for processes with quark-antiquark and gluon-gluon 
collisions by comparing our predictions to the full NNLO corrections for Drell-Yan and Higgs 
production. 

The quantities da^, Po, K, (2, and (3 have all been defined in Section 2.1, and (4 = 7r^/90. 

(2) 

We note that ReF'^. is the real part of the two-loop Ts after dropping all gauge- dependent 
terms. The universal part of 2Rer'^'' — J2i '^/f in quark-antiquark and gluon-gluon collisions 
will be derived in Section 3 from Drell-Yan and Higgs production. B'^^^ will also be derived 
in Section 3. There is also current work on full two- loop evaluations for general processes 
|24j . Also note, again, that the sum over i involves summing over the two incoming partons in 
hadron-hadron collisions, or one parton in lepton-hadron collisions, and that there is no sum 
over j if there are no massless partons in the final state at lowest order. Also note that 6k in 
the formula again indicates which terms vanish in PIM kinematics. Finally, we note that the 
choice = s further reduces the number of terms in the master formula in IPI kinematics; 
however, we keep our expression as general as possible. 

As we will see below, the NNLO master formula passes many rigorous tests. It reproduces 
the exact NNLO soft plus virtual results for Drell-Yan and Higgs production, and deep inelas- 
tic scattering, as well as the NNLO-NNLL results that have been derived already for many 
different processes from threshold resummation studies. Also I have checked explicitly that 
at NNLO the renormalization and factorization scale dependence in the physical cross section 
(after convoluting the partonic cross section with the parton distributions) cancels out for both 
hadron-hadron and lepton-hadron collisions. 



2.3.2 DIS scheme 

In the DIS scheme the NNLO soft and virtual corrections are 

.(2) 



O DIS — " Z2 ^ DIS 



(2.18) 



with 



",(2) 




-,(2) 


^ DIS 




a' \ 


+ 


_ 4 


in- 


+ i 


u 


{T[- 



^'"^:-ET^/.^^2(x,.) 



K 



■ 4 



T^iixth) 



^1. rz^iU^in 



+ { - E 1^. (^) - E Sk f (^) - E Sk f m 



(r^-T2) + E4f}^o(x,.) 

^ In' ( ^\ 



m (^) in - Ti) + 1 in - TO m (^) + E in3 (^) 



-^6k in - T2) In' 



E^f^A'ln 



^ixth) 



(2.19) 
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We note that the T/'s are the DIS quantities while the unprimed Tj's are the MS counterparts. 

Also a' \c'. denotes the cross section in Eq. ()2.17|) after replacing all the Cj, Tj, and R, by 
their DIS counterparts c'j, T/, R'. Our formula reproduces the exact known NNLO corrections 
for Drell-Yan and W'^'y production and deep inelastic scattering in the DIS scheme, as we will 
see in the next section. We will derive the expression for B^^^ in the next section by matching 
to the known NNLO corrections for the Drell-Yan process. We also note that again we have 
checked the scale independence of the physical cross section at NNLO. 



3 Electroweak/Higgs, QCD, and SUSY processes with 
simple color flows 

We now apply our NLO and NNLO master formulas to a variety of processes with simple color 
flow which are of electroweak, QCD, or SUSY origin at lowest order. 



3.1 The Drell-Yan process 

Our first application is the Drell-Yan process, i.e. lepton pair production in hadron-hadron 
collisions, for which the NNLO corrections have been calculated in Refs. jHl IH CS]. Also a 
comparison of the expansion of the resummed cross section in |2IIj with the NNLO corrections 
in [TUj was presented in Ref. The partonic process we discuss is qq ^ V + X, where is a 
vector boson (7, Z, W) which later decays to a lepton pair V —>■ I1I2 with invariant mass 
and X denotes any additional partons in the final state. At threshold s = Q^, where s is the 
center-of-mass energy squared of the incoming quark-antiquark pair. The NNLO corrections 
for the cross section da/dQ"^ in the MS scheme are given in Ref. The soft and virtual 
corrections are given the label A^^^'"^^^ and are given explicitly at NLO by Eq. (B.3) in Ref. 
0, and at NNLO by Eq. (B.8) plus the renormalization scale term in Eq. (B.7) of Ref. 
Here the plus distributions are Di{x) with x = Q'^/s. 

We are able to reproduce these results using our master NLO and NNLO formulas. Evidently 
at lowest order there are no final-state massless partons and the cross section is given in PIM 
kinematics. We note that for this process ReF'^^'* = Cp and we choose the hard scale = Q^. 

At NLO our master formula, Eq. (j2.9j) . gives 

4lv = <,^y^^ {csV^ix) + cMx) + ci5(l - x)} (3.1) 

TT 

with C3 = 4:Cf, 

,._2C.ln(|). cf^-|c, (3.2) 

which reproduces Eq. (B.3) in Ref. j^, and we identify the non-scale 6{1 — x) terms as 
Ti = 2CfC2 - 4:Cf. 

At NNLO our master formula, Eq. ()2.17|) . reproduces the D^i^x), D2{x), Di{x), Do{x), and 
5(1 — x) terms. Note that our results use explicitly the beta function Pq, thus simplifying the 
expression in Eq. (B.8) of pi. We also note that we use the full NNLO soft-plus- virtual result 
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in jH] to derive the two-loop soft anomalous dimension, that appears in the Do{x) term in our 
master formula, for the Drell-Yan process : 

Rer'g), - ) = C,C, (^Ca + yC. - ^) + n.C, (-^G + |) • (3.3) 

Note that this term is universal for all processes with quark-antiquark annihilation. We also 
determine for the Drell-Yan process in the MS scheme the term R that appears in the 6{1 — x) 
terms in our master formula and which is also universal in all qq processes: 

^2/ 59.2 29. 15. 1\ ^ ^ / 3 .2 37. 7. 1535\ 

R,, = ^H"l0^^ + ¥^^"T^' + ^^^'~64j+^^^H"20^^ + y^^ + 4^'~l92"j 

The NNLO soft and virtual corrections for Drell-Yan production have also been calculated 
in the DIS scheme in Ref. 0. Using our NLO formula in the DIS scheme, Eq. ()2.13|) . we find 
c'3 = 2Cf, c'2 = 3Cf/2 - 2CF\n{fil/Q'), c'^ = -i3/2)CF\n{fil/Q^), and T[ = ACfC2 + Cf/2, 
which agress with Eq. (A. 3) in 

Using our NNLO master formula in the DIS scheme, Eq. ()2.18j) . we are also able to rederive 
the NNLO result in Eq. (A. 8) (plus the renormalization scale terms in Eq. (A. 7)) of Ref. j^, 
and thus we identify the B^^^ term in the DIS scheme master formula: 

r.(2) ^2/3 3^ 3 \ ^ ^ / 5^ 4937\ 409 ^ ,^ , 

=Cl{ C2 + -Cs + CfCa — Cs H nfCF ■ (3.5) 

« ^V32 4^ V 2^'^ 864 / 432 ^ ^ ' 

Finally, we also determine for the Drell-Yan process in the DIS scheme the term R' that appears 
in the 5(1 — x) term in our master formula: 

0/ _ r^2 f 43 2 17 . 9 1\ / 77 2 1049 49 215\ 

R,, - ^f[-^C2-j^C2 + -C3 + 3C,--)+CfCa[^--C, + ^G--Cs + j^) 

/Cs 85 19\ 

This is also universal in all qq processes in the DIS scheme, as we will verify below for W'^'y 
production. 



3.2 Standard Model Higgs production 

Our next application is Higgs production in the Standard Model in hadron-hadron collisions, for 
which the full NNLO corrections, using an effective Lagrangian for the Higgs-gluon interaction, 
have been calculated in Refs. El EI- The partonic process we discuss is gg ^ H + X , where 
H is the Higgs boson. At threshold s = Mjj, where s is the center-of-mass energy squared of 
the incoming gluon pair and Mh is the Higgs mass. Evidently at lowest order there are no 
final-state massless partons and the cross section is given in PIM kinematics. Here the plus 
distributions are Di{x) with x = Mfj/s, ReF'^^'' = Ca, and we choose the hard scale = Mfj. 
The soft and virtual corrections for the total cross section Ogg^H are given explicitly at NNLO 
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in Refs. j2l El El IHl 1121 ■ reproduce and generalize those results by keeping the factorization 
and renormalization scales distinct, using the beta function jSo explicitly, and keeping the color 
factors Ca explicit in our results. 

At NLO, our MS scheme master formula gives 



'^SU = ^f,-,H^^ {csVii^) + C2Vo{x) + c,6{l - x)} (3.7) 



IT 



with C3 = 4Ca, 




c, = -2C^ln(^) , c^ = ^ln(^), (3.^ 



which reproduces the results in 00 01111 113 • We also identify Ti = ll/2 + 2(:7^C2- At NNLO, 
our master formula reproduces the D^i^x), D2{x), Di{x), Do{x), and 6{1 — x) terms. We note 
that we use the full NNLO soft-plus-virtual result in these references to derive the two-loop 
soft anomalous dimension, that appears in the Dq{x) term in our master formula, for Higgs 
production: 

Rer'g' - ' = C-c, + ^C. - 4) + 'vc. ( - 4) ■ (3.9) 



This anomalous dimension is universal for all processes with gluon-gluon fusion. 

We also determine for Higgs production in the MS scheme the term R that appears in the 
5(1 — x) terms in our master formula and which is also universal in all gg processes: 

9221 67^ 1089^2 165^ / 1189 5^ 

Raa = + — C2 Cl Cs + 108 C4 + C2 + -Cs • (3.10) 

39 2 20 ^ 4 V 144 3 6 / 

3.3 Deep inelastic scattering 

Our methods can also be applied to the coefficient functions in deep inelastic scattering, 7*13' — > 
q, where the distributions are Di{z). We note that here we have a massless parton (quark) in the 
final state. In the MS scheme, our NLO formula gives: C3 = Cp, C2 = — 3(7^7/4 — Ci7'ln(/i|,/(5^), 
and ci = -{3CF/4:)\n{fil/Q^) - CfC2 - 9Cf/4. The NNLO corrections are then given by 
our NNLO master formula. Our NLO and NNLO corrections agree with the results for the 
coefficient functions in Refs. [5*, 'W\ (see Appendix B in and Appendix A in jH]) and we 
identify the two-loop B'^ that appears in the Dq{x) terms in our master formula: 

^ ^ V32 4^ 2^V V 12^ 864/ V 6^ 432/ ^ ' 

I have checked that in the DIS scheme the NLO and NNLO corrections that do not involve the 
scale vanish, as expected (after all this is the definition of the DIS scheme, that the corrections 
to deep inelastic scattering in that scheme vanish). This involves checking that 2T'^g^^q- 
B'f + (/?o/4)L'g + Bf = which is a further test of the correctness of the expressions for the 
various two-loop quantities that we have derived. 
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3.4 VF^7 production 

We now discuss the cross section s^d?a/{dtdu) for production in pp collisions for which 

the NNLO soft-plus- virtual corrections in the DIS scheme have been presented in Ref. j^Hl- The 
lowest-order partonic process is qq — > and ReV'f = Cp. We define S4 = s + t + 

with s = {j)q+ Pq)^, t = {pq — p-y)^ , u = {j)q — p-y)^ , aud ti = t — iTi^r , Ui = u — m^. The plus 
distributions are Di{si), and we choose = s. 

In the DIS scheme at NLO, our master formula gives 



(1) 



a 



TT 



{C3X>1(S4) + C2X>o(s4) + 4^(54) } 



with Co = 2Cf 



and 



2 



tiUi 



Cl 



3 , 

h In 

2 




Our NLO expansion agrees with Eq. (3.2) in [221 and we identify 



T[ = -Cf In^ 



-Ui 



tiUi 



+ 4:CfC2 + 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



Using our NNLO master formula in the DIS scheme, we are also able to rederive the NNLO 
soft and virtual corrections for that process in Eqs. (Bl) and (B2) of Ref. ^|. Since this is 
the only process previously calculated at NNLO not involving PIM kinematics, and hence the 
ln(— tj/M^) terms are explicit, this re-derivation provides an additional highly non-trivial test 
of our NNLO master formula. 

Finally, we note that we can easily derive the NNLO corrections for this process in the 
MS scheme for which no results have been given in the literature. Now, C3 = 4(7^?, C2 = 
-2C^ln(tiMi/s2)-2Cirln(/i|/s), dj = and Ti = C^. In2(-ti/s) -|- Cj. In^j-wi/s) + 2CK2 - 
ACf- The NNLO corrections in the MS scheme are then given by our master formula. 



3.5 Charged Higgs production 

We now present new results for processes for which no full NNLO calculations have ever been 
done. 

We first consider the production of a charged Higgs boson in the Minimal Supersymmetric 
Standard Model, for which there have been recent NLO calculations [27||2H1- The lowest-order 
partonic process is b{p{) + g{pg) — > H'^{pfj+) + t{pt). Again there are no final-state massless 
partons and in this case we work in IPI kinematics. Here a can stand, for example, for 
Eh+ da/d^pH+- The relevant hard scale we choose is M = mH+ and we define the Mandelstam 
invariants s = {pi + Pg)'^, t = {ph+ — Pi)^, and u = {ph+ — PgY- The threshold variable is 
S2 = s + t + u — rn\j+ — — and the plus distributions are Di{s2). Also, tg = ti = t — 

and ti = Ui = u — mf. The one-loop soft anomalous dimension is ReF'^^^ = Ci?[ln(— -Ui/s) + 
(l/2)ln(s/m2)] + (C^/2)[ln(tiM) + 1]. 
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The NLO soft-plus-virtual corrections in the MS scheme are 



R) 



{C3P1(S2)+C2P0(S2)+C15(S2)} 



with C3 = 2{Cf + Ca), 



C2 = Cb 



In 



H 



smt 



Ca 



In 



m 



H 



, - In — 

tiUi I \ s 



and 



{mir \mir 4 4 



ln(^)+^ln 



(3.16) 



(3.17) 



(3.18) 



The NNLO soft and virtual corrections are then explicitly given by our master formula. 



3.6 VF, Z plus jet production 

Z plus jet production in hadron colliders has been studied at NLO in Refs. [211 EO], and 
a resummed cross section and NNLO-NNLL corrections have been presented in Ref. |19j . 
Here we follow the notation of Ref. [121 ^^d discuss the IPI cross section Eg da/cPQ with Q 
the momentum of the electroweak boson. At lowest order there are two partonic processes, 
(l{Pa) + g{Pb) q{Pc) + V{Q) and q{pa) + q{pb) g{Pc) + V{Q), where V stands for W or Z. 
We note that we have final-state massless partons in both processes. We define the kinematic 
invariants s = {pa +PbY, t = {Pa — QY, and u = {pb — QY- We choose = Q^, and the plus 
distributions are Di{s2) with S2 = s + t + u — Q"^. We discuss the MS corrections for the two 
partonic processes in turn. 



3.6.1 qq^gV 



Here tq = u, tq = t, and ReF^^^ = Cf + {Ca/2) \n{tu/s^) + CA/2 The NLO soft and virtual 
corrections are 



^Sl.v = ^!,-.,v^^ {csV^is^) + C,V,{S,) + C15(S2)} 



TT 



with cs = ACf — Ca, 



C2 = -9^ - 2Cpln ^ - {2Cf - Ca) In 



tu 



and 



di = CF 

These are in agreement with the NLO result in [^ . 



3 , tu 



.(^).§.(^). 



(3.19) 



(3.20) 



(3.21) 
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3.6.2 qg qV 



Here tg = u, tg = t, and ReT'^s^ = CpH-u/s) + Cp + (Ca/2) ln{t/u) + Ca/2 HH]. The NLO 
soft and virtual corrections are 



.(1) 



^qg^qV - ^qg-^qV 



n 



{C3P1(S2)+C2P0(S2)+C15(S2)} 



with c3 = Cf + 2C/ 



(3.22) 



C2 



-Cf - {Cf + Ca) In 



-Ca In 



tu 



(3.23) 



and 



A_^CF + CF\n (^]+CA\n(=^' 



(3.24) 



They are in agreement with the NLO result in [HU] . 

For both partonic subprocesses the NNLO corrections are derived from our master for- 
mula and are in agreement with the NNLO-NNLL results in Ref. [19j (apart from a term 
(j'^[Q!f(/i|.)/7r^](/3o/4)c3ln((5^/s)Pi that was missing in that reference). 



3.7 Direct photon production 

Direct photon production is often recognised as a process that can aid determinations of the 
gluon distribution. The NLO cross section for direct photon production has been given in 
Refs. |nillS2]- At lowest order, the parton-parton scattering subprocesses are q{pa) + g{pb) 
liP-y) + i{Pj) and q{pa) + giPb) liP'y) + 9{.Pj)i so there are final-state massless partons in 
both subprocesses. We define the Mandelstam invariants s = {j)a+ PbY, i = ij>a — P-yYi and 
u = {pb—p-yY, which satisfy S4 = s+t+u = at threshold. Here we choose = Pt = tu/s, and 
we work in IPI kinematics in the MS scheme with the cross section d^a/(fp^. The threshold 
logarithms -D/(s4) have been resummed and NNLO-NNLL corrections have been presented in 
Ref. (33j. 



3.7.1 qq ^ 75f 

We start with the process qq 7(7 for which tg = u,tq = t, and ReF'^^^ = Cf+{Ca/2) ln(tn/s^) + 
Ca/^,. The NLO soft plus virtual corrections are 

41^, = af^^^,^^ {csV.is,) + C2Vo{s,) + ci5(s4)} (3.25) 
with C3 = ACf — Ca, 

C2 = -§-2C.ln(^) , (3.26) 
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and ci = Ci + Ti with 



•-1 



2 Is, 



(3.27) 



The term Ti is given in Eq. (3.11) of Ref. [22] (where it is called c'^''). The NNLO soft plus 
virtual corrections are then explicitly given by our master formula. We note that the terms 
-03(54), -02(54), and -Di(s4) were already derived from a resummation study in Section III of 
Ref. [22] and are in agreement with our formula. 



3.7.2 qg -fq 

We continue with the process qg — > 7g for which tq = u, tg = t, and Rer'^^"* = Cp^Jii—u/ s) + 
Cp + (C*a/2) \n{t/u) + Ca/'^- The NLO soft plus virtual corrections are 



cr 



(1) 



a 



B 



R) 



{oiVi{si) + 02^0(54) + ci5{sa)} 



with C3 = + 2Ca, 



and Ci = Ci + Ti with 



/3o 3 



CfIu 



-t 



C.ln 



— M 



In 



+ -ln 



(3.28) 



(3.29) 



(3.30) 



The term Ti is given in Eq. (3.8) of Ref. |22] (where it is called c'^^). The NNLO soft plus 
virtual corrections are then explicitly given by our master formula. Again, we note that the 
terms .03(54), .02(54), and 0*1(54) were already derived from a resummation study in Section 
III of Ref. [22] and are in agreement with our formula. 



3.8 DIS heavy quark production 

The NNLO corrections to heavy quark production in deep inelastic scattering may be needed, 
together with various resummations |Ml l35j. in explaining the discrepancy between NLO theory 
PU] and experiment for bottom quark production. Here the lowest-order partonic process is 
1*9 ~^ QQi so there are no final-state massless partons, and we are working in IPI kinematics in 
the MS scheme with the cross section (Pa/{dtidui). The soft anomalous dimension is ReF'^^'* = 
(Ca/2 - CF){ReLf3 + 1) + (Ca/2) \n{tiUi/{mh)), where 5 = (p^. +Pgy, h = {Pg - PqY - 
and Ml = (p^« —pgY — m^, with m the heavy quark mass. The singular distributions are Di{s2) 
with 52 = 5 + ti + Ml, and we use M = m. 
The NLO corrections are 

= {c,V,{s,) + c,V,{s,) + ci5(52)} (3.31) 
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with C3 = 2Ca, 



C2 = -2CF{ReLf^ + 1) + Ca 



ReL/3 + In 



In 



and 



In 



+ -ln 



(3.32) 



(3.33) 



The NNLO soft and virtual corrections are then given by our master formula. We note that 
the 1^3(52) and 02(32) terms were derived in Ref. jSl] (with fip = fi^i) and are in agreement 
with our results. 



3.9 FCNC single-top production 

The last process with simple color flow that we consider is single-top production mediated 
by flavor-changing neutral currents (FCNC). The QCD corrections to the IPI cross section 
da/{dtdu) in the MS scheme for FCNC single-top-quark production in ep collisions at the 
HERA collider were studied at NLO in the eikonal approximation in Ref. Here the partonic 
process is eu —>■ et, so there are no final-state massless partons, and ReF'^^'* = CplnKm'f — 
t)/(-\/smt)], with s = {pe + PuY, t = {Pt — PuY, and u = {pt — PeY- The singular distributions 
are Di{s2) with S2 = s + t + u — — 2ml, and we choose M = rrit. 
The NLO corrections are 



;^(l) B (^sijJ^r ^ , 

71 



S2) + C2'Do{s2) + Ci6{s2)} 



(3.34) 



with C3 = 2Cf, 



C2 



-1 - 21n 



-u 



+ 2 In 




(3.35) 



and 





(3.36) 



These agree with the results in |3ZI. The NNLO corrections are then given by our master 
formula. 



4 NNLO master formula and applications - complex color 
flow 

4.1 NLO and NNLO master formulas for complex color flow 

When the lowest-order cross section involves already a complex color flow that is expressed 
in terms of non-trivial color matrices, then the master formulas given in Section 2 have to be 
extended. Now not all the NLO corrections are proportional to the Born term; only the leading 
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logarithms and terms involving the scale are. Therefore at NLO our master formula for soft 
and virtual corrections in the MS scheme is extended for the case of complex color flow as 



cr 



(1) 



R 



a 



R) 



TT 



{c^Vi{xth) + C2VQ{xth) + ci5{xth)} + 



a 



TT 



where C3 is defined as before in Eq. ()2.1()|1 . and C2 is defined by 



(4.1) 



C2 



E 



Cf^ + 2Cf^ 5k In 



E 



B'f + Cf^ + Cj^ 5^ In ( — 



M2 



(4.2) 

i.e. is the same as in the simple color flow case except without the term 2Rer'5''. The function 
is pro cess- dependent and depends on the color structure of the hard-scattering, for which 
specific examples are given in the next subsections. It is defined by 



A' = Tr (i7Wr«t^(o) + ff(o)5(o)r/a)^ _ 



(4.3) 



Note that we use the expansions for the hard and soft matrices: H = H^^'> + {ag°'^~^^ / 7r)H^^^ + 
(af »+V7r2)i7(2) + . . . and ^ = + (a./vr)^^!) + (a./vr)^^^^) + ■ ■ The Born term is then 
given by = Tr[if(o)5(o)]_ 

With respect to the 6{xth) terms, we split them into a term ci = Ci + Ti, with Ci defined in 
Eq. ()2.12j) as before, that is proportional to the Born cross section, and a term Tf that is not. 
Tf is also pro cess- dependent, is formally defined by 



Tr fi7(^)^(°) + H^'^S^'A + A^6k In ( —] + a^"^^^^) 



TT 



C3 



+ T2 5k In 



(4.4) 



and can also be derived by matching to a full NLO calculation. 

In the DIS scheme the corresponding terms C3, Cg, and c'l are given as in Section 2.2.2. 

At NNLO the master formula for soft and virtual corrections in the MS scheme is extended 
for the case of complex color flow as 



a 
+ 
+ 



(2) 



"simple ' 



a 



(/^l) j 3 



77^ 



-C3A'V2{Xth) + 



'202-^] A' + c^T^ + F 



;c.-C.C3.f,n(f 



+ C2 - 



T^ + F6k In 



M2 



(C3C3 - (202)^^= + c, + ^\n r^ T^ + - 6K\n' — 



M' 



A 



^SKT^lnl — ] +^SKA'\n^{ — \ +G5k In ( — 



M2 



M2 



-C2 



T'liXth) 

G V,{xth) 



(4.5) 
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fry's. ^ ^ 

Here (Tgij^pig denotes the expression in Eq. ()2.1fi|l after using the new C2 of Eq. ()4.2j) everywhere 
in that expression and deleting all T's, and R, from that expression. Also, we have used 



Tr 



(4.6) 



and 



G 



Tr 



^(i)r/(i)t^(o) ^ ^(D^Wr'W + ijWr'Wt^d) 



^(0)^(1)^/(1) 



^ ^(0)^/(2) t^(0)^^(0)^(0)r/ 



i?^ = Tr (i/(2)^(o) + Hio)si2) ^ jjiDgii 



2 1 



(4.7) 



(4.8) 



I have checked explicitly that the renormalization and factorization scale dependence cancels 
in the physical cross section at NNLO. 

In the DIS scheme the NNLO corrections are given by Eq. (j2.18|) . after replacing the term 

{aKfij^) / 'n''^)a' \c'. by the cross section in Eq. ()4.5|) . having also replaced in Eq. fl4.5|) all the 
Cj, Tj, G, and i?*^, by their DIS counterparts c^, T/, G', and W^' . 

We now apply the NLO and NNLO master formulas to a variety of processes with complex 
color flow. 



4.2 Heavy quark hadroproduction 

The production cross sections of top, bottom, and charm quarks in hadron colliders can be 
considerably enhanced near threshold. The NLL resummation of threshold corrections was 
derived in Ref. ^3] and the NNLO-NNLL corrections for heavy quark total cross sections and 
differential distributions were calculated in Refs. fSl IHH] in both IPI and PIM kinematics. 
There are two partonic channels involved at lowest order, qq — >■ QQ and gg QQ. We choose 
the hard scale M = m, with m the heavy quark mass. The IPI cross section is d?a /{dti dui), 
with ti = t — m^, ui = u — m^, and the singular distributions are Vi{s4) with S4 = s + ti + ui. 
The PIM cross section is s d'^a/ldM^Qd cos 6), with Mqq the heavy quark pair mass and 6 the 
scattering angle in the partonic center-of-mass frame, and the singular distributions are 'Dilz) 
with z = Mqq/s. The one-loop soft anomalous dimension matrices P^'' were calculated for 
both partonic processes in Ref. ^3]. Explicit results for the matrices H^^^ and S^^^ can be 
found in Refs. fSl IHHI for both partonic channels. 

4.2.1 qq ^ QQ channel 

We begin with the quark-antiquark annihilation channel. We note that P'^, H, and S are 2x2 
matrices. However, the triviality of the hard matrix H^^^ (only one non-zero element) leads to 
almost simple-color-flow-like expressions. 

The NLO soft plus virtual corrections in IPI kinematics in the MS scheme are given by Eq. 
(jUD with C3 = AGf, 

c, = -2Gf - 2Gf In - 2Gf In (^^ j , (4.9) 
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^1 



In 



ln(^)+|ln(^^ 



(4.10) 



and A" = {a^ /a^) 2ReT'^sl2^ and are in agreement with the NLO results in jSH]. Here the 

real part of the one-loop 22 element of the soft anomalous dimension matrix is Rer'^^22 — 
Cir[41n(Mi/ti)-ReL^] + (CA/2)[-31n(Mi/ti)-ln(m2s/(tiMi))+ReL^] with = (l-2mVs)//3- 

[ln((l — /?)/(l + /?)) + vri] and /? = y^l — im'^/s. Exphcit expressions for Ti, Tf can be extracted 
from Ref. [39 . 

The corresponding NLO MS corrections in PIM kinematics are given by similar expressions 
after striking out the ln{tiUi/m'^) terms from the IPI C2 and c^, as per our NLO master formula, 
and using the relevant PIM Born term and Ti, |38j . 

The NNLO soft and virtual corrections are then given by our master formula for complex 
color flows for either kinematics and are in agreement with the NNLO-NNLL results in [H^ I38j . 
We note that the term F in Eq. ()4.6|1 has a relatively simple form, F = ((j'^/a;g)[4(Rer'g22)^ + 



4r' 



(1) 



s.is^'s^zi] "w-ith 



iCF/CA)Hui/ti) and r^sl, = 2\n{ui/h). 



In the DIS scheme, we have in IPI kinematics Cg = 2Cf, 



Co = C i? m 

^ 2 




(4.11) 



Cl 



2 



In^ 



-Ui 



--cm 



-2CfC2 + \Cf. (4.12) 

The corresponding corrections in PIM kinematics are again given by similar expressions after 
striking out the ln^(— ti/m^), ln^(— Wi/m^), and ln(tiMi/m^) terms from the IPI c'2 and c'^^, and 
using the relevant PIM Born term and T^' [SB]- The NNLO soft and virtual corrections in the 
DIS scheme are then given by our master formula for either kinematics. 

4.2.2 gg QQ channel 

We continue with the gg QQ channel whose color structure is considerably more complex. 
Here P'^, H, and S are 3x3 matrices. The NLO soft plus virtual corrections are given in IPI 
kinematics in the MS scheme by Eq. (I4.1|l with C3 = 4Ca, 



C2 = -2Ca - 2C^ln ('^) - 2C^ln 



(4.13) 



and 
A" 



CaIii 



T^K,gBQED{Nl-l){N, 



tiUi 



2 



tl + ul 



c 



+ (ReL^ + l) + -^ + -^ln 



+ ^(C.-C.)(ReL, + l)-ln(^)+ ^ 




(4.14) 



ni?s 



(4.15) 



20 



and are in agreement with the NLO results in IIU]. Exphcit expressions for Ti, Tf can be 
extracted from Ref. [lU]- As for the qq channel, the corresponding NLO corrections in PIM 
kinematics are obtained after striking out the ln(iiMi/m^) terms from the IPI C2 and Ci (note 
that is not affected by the kinematics choice) and using the relevant PIM Born term and 

Ti, inni- 

The term F can be explicitly calculated using Eq. ()4.6|) . The NNLO soft and virtual 
corrections for either kinematics are then given by our master formula for complex color flow 
and are in agreement with the NNLO-NNLL results in [TH| I38j. 



4.3 Jet production 

Threshold resummation for jet production has been studied in Refs. fHl E]- There are many 
partonic subprocesses for which explicit results for the NNLO corrections at NLL accuracy 
were given in |4^. Here we are able to extend those results by employing our master formula. 
We note that the color structure for these processes is quite complex, and the one-loop soft 
anomalous dimension matrices fHl ^] along with the lowest-order hard and soft matrices can 
be found in Ref. [H] . Recently, the complete two- loop virtual corrections have been calculated 

miaiiaiisi. 

Here we discuss the single-jet inclusive cross section Ej d^a/d^pj in the MS scheme. The 
NLO soft and virtual corrections with = p\ = tu/ s and S4 = s + t + m can be written for 
each subprocess as 

^(1) ^ |C3P,(,^) + ^2^0(34) + ci5{s,)} + [A'^V,{s,)+Tt5{s,)] . (4.16) 

TT TT 

The expressions are symmetric in t and u except for the qg qg channel. Full NLO expressions 
have been given in [Hj. We discuss next the individual partonic processes in jet production. 



4.3.1 qq qq and qq — > qq 

There are several processes involving distinct or identical quarks and antiquarks: qjqj qjqj, 
IjQj ~^ QkQk, QjQk QjQk, QjQj Ijlj, QjQk QjQk, and the corresponding ones with antiquarks. 
The NLO soft and virtual corrections are given by Eq. ()4.16|) with C3 = 2Cf, 

c, = -2CF\n(^)-^CF, (4.17) 




The expressions for cr^ and depend on the specific process. The expressions for can be 
found in Appendix A (for qq qq processes) and Appendix B (for qq qq processes) of Ref. 
|41j . The expressions for A'^ can be easily derived from Eq. (|4.3|) or by comparing Eq. (j4.16p 
with the expressions in Appendix A or Appendix B of [33] for the NLO corrections for the 
various subprocesses. 
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4.3.2 qq gg and gg qq 

The NLO soft and virtual corrections are given by Eq. ()4.16|) with C3 = ACp — 2Ca, 



C2 = -2Cf In 



2 



2Cf 



2Ca - 2Ca In 



for the process qq — >• (75^, and C3 = 4Cyi — 2Cf, 



c, = -2CA\n 



]^Cf-2Ca-2Cf\^ 



(4.19) 
(4.20) 

(4.21) 
(4.22) 



for the process gg qq. The expressions for depend on the specific process and can be 
found in Appendix C of Ref. [H]. The expressions for can be easily derived from Eq. ()4.3|) 
or by comparing Eq. ()4.16p with the expressions in Appendix C of 01] for the NLO corrections 
for these subpro cesses. 



4.3.3 qg qg 

The NLO soft and virtual corrections, using tg = u and tg = t, are given by Eq. ()4.16|) with 

C3 = Cf + Ca, 

C2 = -{Cf + Ca) In f 41 - ^^f -2Ca-^- 2C,.ln - 2C^ln , (4.23) 



and 



Cpln{-l^+CAln (^-^ 



4^^ + T 



ln(^)+|ln(^^ 



(4.24) 



The expression for can be found in Appendix D of Ref. [H] . The expression for A'^ can be 
easily derived from Eq. ()4.3j) or by comparing Eq. ()4.1fj|l with the expression in Appendix D 
of 03] for the NLO corrections for this subprocess. 



4.3.4 gg gg 

The NLO soft and virtual corrections are given by Eq. ()4.1()j) with C3 = 2Ca, 

c, = -2CAln (^)-y-4C^A, 



(4.25) 



and 



(4.26) 
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The expression for can be found in Appendix E of Ref. |3I]. The expression for A'^ can be 
easily derived from Eq. ()4.3|) or by comparing Eq. ()4.16p with the expression in Appendix E 
of |H] for the NLO corrections for this subprocess. 

The term F in Eq. (j4.6p can be easily derived for each subprocess using the matrices 
H^^\ 5'(o)^ and T^g^ as given in Ref. [H]. For all subprocesses the NNLO soft and virtual 
corrections are given by our master formula for complex color flows and are in agreement with 
the NNLO-NLL expressions in [H]. 



4.4 Squark and gluino production 



Our last application is the production of squarks and gluinos in hadron colliders. The complete 
NLO corrections for squark and gluino production in the Minimal Supersymmetric Standard 
Model have been given in [33]. There are several partonic subprocesses involved. 

We begin with squark production. For the channel qq —>■ qq the NLO terms C3, C2, and 
Ci in the MS scheme are the same as for the qq QQ channel in heavy quark production in 
Section 4.2.1 in either IPI or PIM kinematics, with M = m the mass of the squark. The same 
holds for the channel qq qq, and analogous results hold in the DIS scheme. For the channel 
99 ~^ QQ ill the MS scheme, C3, C2, and Cj* are the same as for the gg QQ channel in heavy 
quark production in Section 4.2.2 in either kinematics. 

We continue with gluino production. Here we choose M = m, with m the mass of the 
gluino. For the process qq gg, C3, C2, and are the same as for the process qq qq in 
either IPI or PIM kinematics and either factorization scheme; for the process gg gg in the 
MS scheme they are the same as for gg —>■ qq in either kinematics. 

We finally have the process of squark-gluino production, qg — > qg. We choose M = m, with 
m the mass of the squark or the mass of the gluino. In IPI kinematics in the MS scheme we 
have C3 = 2{Cf + Ca), 



C2 



-Cf-Ca-2Cf\^ 



-Ml 



{Cf + Ca) ln(i^' 



(4.27) 



and 



C^ln 



-Ui 



CAln 



3^ 



(4.28) 



In PIM kinematics the expressions for C2 and are obtained after striking out the ln(— ti/m^) 
and ln(— Ui/m^) terms from the IPI C2 and Ci. 

For all these processes, A"^, Ti, can be read off the full NLO calculation. Also A^ and F 
can be calculated explicitly once the lowest-order H, S, and Ts matrices have been constructed. 
The NNLO soft and virtual corrections are then given by our master formula. 



5 Conclusions and outlook to higher orders 

In this paper, I presented a unified approach to calculating the NNLO soft and virtual QCD 
corrections for any process in hadron-hadron and lepton-hadron collisions in either the MS or 
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DIS schemes and in either IPI or PIM kinematics. The master formulas given in the paper are 
based on a unified threshold resummation formalism and they allow explicit calculations for 
any process, with either simple or complex color fiows, keeping in general the factorization and 
renormalization scales separate and the beta function and color factors explicit. I verified that 
the scale-dependence of the physical cross section cancels out at NNLO for any process. Detailed 
results, illustrating the use of the master formulas, were given for various electroweak, Higgs, 
QCD, and SUSY processes in various factorization schemes, kinematics, and colliders. As tests 
of the master formulas, I reproduced the previously known NNLO corrections for Drell-Yan 
and Higgs production, deep inelastic scattering, and W'^'j production, thus also determining 
a number of two-loop anomalous dimensions and other quantities which are needed in NNLL 
resummations. Furthermore, I presented new results for several other processes in the Standard 
Model and beyond. 

The NNLO corrections increase theoretical accuracy and diminish the dependence on the 
factorization and renormalization scales, and thus are essential in further testing QCD and 
particularly in searching for the Higgs boson and supersymmetric particles as well as other 
processes, such as fiavor- changing neutral currents, that signal new physics beyond the Standard 
Model. 

The unified approach employed in this paper can be extended to higher orders. For a sketch 
of how this extension may be carried out through next-to-next-to-next-to-next-to-leading order 
in the specific context of heavy quark hadroproduction see Ref. jlH]. A complete calcula- 
tion of next-to-next-to-next-to-leading and even higher-order soft and virtual corrections is a 
formidable task and unlikely to be necessary at least in the forseeable future. The accuracy 
attainable at NNLO should be sufficient for the high-energy colliders of our era. The unified 
master formulas for the NNLO soft and virtual corrections for any process, which are an im- 
portant component of the full NNLO calculation, should serve as a milestone in the push for 
ever-increasing theoretical accuracy and understanding of high-energy processes. 
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